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Diagnostics of irradiated gas in galaxy nuclei 

I: A Far-ultraviolet and X-ray dominated region code 
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Abstract. We present a far-ultraviolet (PDR) and an X-ray dominated region (XDR) code. We include and discuss thermal and 
chemical processes that pertain to irradiated gas. An elaborate chemical network is used and a careful treatment of PAHs and 
H2 formation, destruction and excitation is included. For both codes we calculate four depth-dependent models for different 
densities and radiation fields, relevant to conditions in starburst galaxies and active galactic nuclei. A detailed comparison 
between PDR and XDR physics is made for total gas column densities between ~ 10 20 and ~ 10 25 cm~ 2 . We show cumulative 
line intensities for a number of fine-structure lines (e.g., [CII], [OI], [CI], [Sill], [Fell]), as well as cumulative column densities 
and column density ratios for a number of species (e.g., CO/H2, CO/C, HCO + /HCN, HNC/HCN). The comparison between the 
results for the PDRs and XDRs shows that column density ratios are almost constant up to JVh = 10 22 cm -2 for XDRs, unlike 
those in PDRs. For example, CO/C in PDRs changes over four orders of magnitude from the edge to JVh = 10 22 cm -2 . The 
CO/C and CO/H2 ratios are lower in XDRs at low column densities and rise at Nb. > 10 23 cm -2 . At most column densities 
A^h > 10 21 ' 5 cm -2 , HNC/HCN ratios are lower in XDRs too, but they show a more moderate increase at higher A^h. 



1. Introduction 

Gas clouds in the inner kpc of many galaxies are exposed to 
intense radiation, which can originate from an active galac- 
tic nucleus (AGN), starburst regions or both. O and B stars 
dominate the radiation from starbursts, which is mostly in the 
far-ultraviolet (6.0 < E < 13.6 eV), turning cloud surfaces 
into Photon Dominated Regions (PDRs, Tielens & Hollenbach 
119851) . Hard X-rays (E > 1 keV) from black hole environ- 
ments (AGN) penetrate deep into cloud volumes creating X-ray 
Dominated Regions (XDRs, Maloney et al. ll996l) . For each X- 
ray energy there is a characteristic depth where photon absorp- 
tion occurs. So for different spectral shapes, one has different 
thermal and chemical structures through the cloud. Although 
one source can dominate over the other energetically, (e.g., an 
AGN in NGC 1068 or a starburst in NGC 253), the very dif- 
ferent physics (surface vs. volume) require that both should be 
considered simultaneously in every galaxy. 

In PDRs and XDRs, the chemical structure and ther- 
mal balance are completely determined by the radiation field. 
Therefore, PDRs and XDRs are direct manifestations of the 
energy balance of interstellar gas and their study allows one 
to determine how the ISM s urvives the presence of star s and 
AGN (Tielens & Hollenbach EH1 Boland & de Jon g 1 19821 
van Dishoe ck & Black 1 19881 L e Bourlot et al. 1 19931 Wolfire 
et al. ll993i Spaans et al. 1 19941 Sternberg & Dalgarno fl995l 
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Stoerrzer et a ll 19981: Spaans 1 1996t Bertoldi & Dra ine Il996t 
Maloney et al. 1 19961 Lee et aL[l996; Kaufman et al. 1 19991 Le 
Petit et al. 2002 and references therein). 

PDRs and XDRs have become increasingly important as 
diagnostic tools of astrophysical environments with the advent 
of infrared and (sub-)millimetre telescopes. PDRs emit fine- 
structure lines of [CI] 609, [CII] 158 and [OI] 63 [im; rota- 
tional lines of CO; ro-vibrational and pure rotational lines of 
H2; many H2O lines as well as many broad mid-IR features 
associated with Polycylic Aromatic Hydrocarbons (PAHs). In 
PDRs, the bulk of H 2 is converted into atomic hydrogen at the 
edge and CO to neutral carbon into ionised carbon. XDRs emit 
brightly in the [OI] 63, [CII] 158, [Sill] 35, and the [Fell] 1.26, 
1.64 /im lines as well as the 2 /im ro-vibrational H2 transitions. 
The abundance of neutral carbon in XDRs is elevated compared 
to that in PDRs and the chemical transition s from H to H2 and 
C+ to C to CO are smoother (M alonev et al J 19961) . 

In this paper, we compare a far-ultraviolet and X-ray dom- 
inated region code. For the PDR and XDR, we discuss the 
cooling, heating and chemical processes. Then we show four 
models with different radiation fields and densities, for a semi- 
infinite slab geometry and irradiation from one side without ge- 
ometrical dilution. We conclude with a comparison between the 
column densities, integrated line fluxes and abundance ratios. 
We will apply this tool to the centres of nearby active galax- 
ies in subsequent papers. We would like to point out that these 
codes can be used over a broad range of physical situations and 
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scales, e.g., young stellar objects, planetary nebulae or gas out- 
flow in galaxy clusters. 

2. The Photon Dominated Region model 

The global properties of PDRs are determined by a number of 
physical processes: 

i. Heating through photo-electric emission by dust grains and 
PAHs (c.f. Bakes & Tielens ll994l Weingartner & Draine 
1200 lb . 

ii Heating by FUV pumping of H 2 , follo wed b y collisional 
de-excitation (c.f. Hollenbach & McKeefl979i) . 

Hi. Heating by cosmic rays (c.f. Field ll969h . 

zv. Fine-structure line cooling of [CI] 609, [CII1 158, [OI] 146 
and 6 3 /zm (c.f. Tielens & Hollenbach 1985; Spaans et al. 
I1994I) . 

v. Molecular line cooling by warm molecular gas containing 
CO, H ? , H ? Q , OH and CH (c.f. Neufeld et al l 199.4 Spaans 

vi. Ion-molecule reactions driven by the ionisation degree of 
~ 10~ 4 maintained by the ionisation of carbon in t he FUV 
(c.f. Black & Dalgarno fl977t van Dishoeck & Black[l986). 

vii. The ionisation balance of atomic gas under the influence 
of photo-ionisation reactions driven by FUV photons and 
counteracting recombination and charge transfer reactions 
with metals and particu larly PAHs (c.f. Lepp & Dalgarno 
Il988t Bakes & Tielens ll994l) . 

As one moves into a PDR the extinction along the line of 
sight increases and the impinging radiation field is attenuated. 
Consequently, there are two zones over which the chemical 
composition of the PDR changes in a fundamental way. The 
first fundamental change occurs at the very edge of the PDR 
as atomic hydrogen is converted into H 2 because the Lyman 
and Werner electronic bands that lead to dissociation of the H 2 
molecule in the FUV become optically thick (so-called self- 
shielding). Deeper into the PDR, at about 3 mag of extinction, 
ionised carbon is quickly converted into neutral form as the 
FUV flux decreases due to dust absorption. C is subsequently 
transformed into CO, since the FUV field is reduced by grain 
opacity, H 2 shielding and some CO self-shielding. 

The first few magnitudes of extinction of the PDR are usu- 
ally referred to as the radical region since many carbon hy- 
drides and their ions, e.g., CH, CH+, CN, HCN, HCO+ (and 
also CO + ), reach their peak abundance there, caused by the 
presence of both C + and H 2 and the high (~ 10 2 - 10 3 K) 
temperatures. Ion-molecule reactions take place that lead to 
the formation of a large number of different molecular species. 
Many of the atoms and molecules in (the radical region of) a 
PDR are collisionally excited at the ambient densities and tem- 
peratures, and emit brightly in the mid-IR, FIR, millimetre and 
sub-millimetre. 

The global characteristics of any PDR are defined by a few 
key parameters: 

;. The strength of the impinging radiation field, Go or Ijjv, m 
units of the Habin gil l 1969b or lDrainel(ll978b radiation field, 
respectively, determines the total available radiative flux at 
the edge of the PDR. 



ii. The temperature and the ambient hydrogen density, nn = 
n(H) + 2n(H 2 ), sets to a large extent the pace of the chem- 
ical reactions and the excitation rates of the coolants. 

Hi. The metallicity Z, in units of the solar value Z@, constrains 
the total abundances possible for carbon- and oxygen- 
bearing species and hence influences the chemical and ther- 
mal structure. 

iv. The spectral shape of the impinging radiation field, param- 
eterised by the colour temperature T e g for black bodies or 
the frequency slope for power laws, fixes the distribution of 
photon flux over energy. 

The details about the thermal and chemical processes we use in 
the code are discussed in the appendices. In the rest of the pa- 
per we use Go, the Habing flux, as the normalisation in which 
we express the incident FUV radiation field, where Go = 1 
corresponds to a flux of 1.6 x 1CP 3 erg cm -2 s _1 . 

3. PDR test models 

In this section, we discuss the results for four PDR models in 
which we have varied the radiation field Go and the density 
nu. The models are for a semi-infinite slab geometry, but the 
code also allows for two-sided slab geometries. The adopted 
model parameters are listed in Tabled Models 2 and 4 will be 
shown in a paper by Rollig et al (in prep.), where they are used 
to compare 12 different PDR codes that are commonly used. 
The parameters are listed in Tabled These values are typical 
for the high density, strong radiation field conditions we want 
to investigate in, e.g., a starburst. 



Table 1. Adopted model parameters 



Model 


Go 


Ffuv 
[erg cm -2 s _1 ] 


n H 
[cm -3 ] 


1 


10 3 


1.6 


10 3 


2 


10 s 


160 


10 3 


3 


10 3 


1.6 


1Q 5.5 


4 


10 5 


160 


1Q 5.5 


Sva (km s^ 1 ) 

s d 


2.7 
1.0 







The fixed gas-phase and total abundances we use are given 
in Ta ble |2] The total abundance s are t he average values of 
lAsolund et ail J2004 and Jenkins] J2004 . To calculate the gas- 
phase ab undances, we use the depletion factors calculated by 
I.TenkinsN2004l) . 

3. 1. Heating and cooling 

For both radiation fields and densities, the dominant heat- 
ing source to a column density N-r ~ 10 22 cnr 2 is photo- 
electric emission from grains. In the moderately low density, 
low radiation-field Model 1, viscous heating is about equally 
important in this range. For Model 2 where the radiation field 



Rowin Meijerink & Marco Spaans: Diagnostics of irradiated gas in galaxy nuclei 



3 



Table 2. Abundances 



Species 


A (gas) 




Oj 


J\i (total ) 


He a ' c 


8.5 x 1(T 


-2 


1.0 


8.5 x 10" 2 


C 


1.4 x 10" 


■4 


0.6 


2.5 x 10" 4 


N 


5.2 x 10" 




0.7 


7.2 x 10~ 5 




OA w ln - 
0.4 X 1U 


-4 


u. / 


4. < X IU 


q; 


L.i X 1U 


-6 


n n^ 


0.4 X IU 


s 


6.9 x 1(T 


-6 


0.5 


1.4 x 10 -5 


CI 


5.4 x 10" 


-8 


0.2 


2.4 x 10~ 7 


Fc 


2.0 x 10" 


• 7 


0.007 


2.8 x 10" 5 


r 


o.y x iu 


■8 


U. 1 


z.y x iu 


Kin 

in a 


o.y x iu 


■7 


U,4 


1.0 X IU 


Mg 


Z. / X iu 


-7 


U.Uo 


o.4 X IU 


JNe 








o.y X IU 


Al 








2.3 x 10 -6 


Ar 








1.5 x 10~ 6 


Ca 








2.0 x 10" 6 


Cr 








4.4 x 10" 7 


Ni 








1.7 x 10" 6 



Notes to Tabled 

a : Present in both PDR and XDR chemical network 
b : Present in PDR chemical network 
c : Used to calculate <r pa for XDR 



is increased, it contributes somewhat more than 10 percent. For 
the low radiation-field, high density Model 3, carbon ionisa- 
tion is the second important heating source. In Model 4, where 
the radiation field is increased compared to Model 3, this is H2 
pumping. At high column densities (iVn > 22.5 cm -2 ), [OI] 
63 [im absorption and gas-grain heating are important. For the 
low density PDRs Model 1 and 2, only [OI] 63 [im dominates. 
When the density is increased in Models 3 and 4, gas-grain 
heating is equally important if not dominant. Other heating 
processes contribute less than 10 percent, but are sometimes 
important in determining the thermal balance. 

In all models [OI] 63 /im cooling dominates to Ah = 
10 21 - 5 cm" 2 . In the low density PDRs, [CII] 158 fjm cool- 
ing contributes more than ten percent of the cooling in this 
range, where at high densities, gas-grain cooling is the second 
most important coolant. In the high density, high radiation-field 
Model 4, this contribution can be almost forty percent. Deeper 
into the cloud, [CI] 610 /im and CO line cooling become im- 
portant. H2 line cooling can contribute up to 10 percent to the 
total cooling rate at some point, but is always a minor coolant. 

3.2. Chemical and thermal structure 

The H — > H2 and C + — > C — > CO transitions are quite sharp. 
Their actual location greatly varies, since this is strongly de- 
pendent on density and radiation field. Exposed to stronger ra- 
diation fields, the transitions occur deeper into the cloud, since 
the photo-dissociation rates are larger. At higher densities, the 
transitions occur closer to the surface of the cloud, since the 
recombination rates scale as n 2 . For the same reason, the H + 



and + fractional abundances are systematically higher in the 
low density models. SiO and CS are more abundant and formed 
closer to the surface in the high density models, which is also 
the case for HCO+, HCN, HNC and C 2 H. 

The edge temperatures (see Fig.[6j are affected most by the 
strength of the radiation field when the density is largest. At a 
density of njj = 10 5 5 cm" 3 , the difference is a factor of thirty 
for an increase from Go = 10 3 to Go = 10 5 . In the low den- 
sity case this is only a factor of two. Because of optical depth 
effects, CO cooling is less effective at high column densities. 
For this reason, temperatures rise again at Ah ~ 10 22 cm -2 in 
the low density models. 

4. The X-ray Dominated Region model 

Unlike PDRs, XDRs are mostly heated by direct photo- 
ionisation of the gas, which produces fast electrons that lose 
energy through collisions with other electrons, as well as H 
and H2. These fast electrons collisionally excite H and H2, 
which subsequently emit Lyman a and Lyman-Werner band 
photons, respectively. These photons in turn are capable of ion- 
ising atoms such as C and Si or ionise and dissociate molecules 
such as H2 and CO. 

Compared to PDRs, t he fol lowing processes play a role in 
XDRs (c.f. Maloney et al. ll996h . in part because of the produc- 
tion of UV photons as described above: 

2. Photo-ionisation heating (i.e., Coulomb heating with ther- 
mal electrons) dominates by a large factor over the heating 
through photo-electric emission by dust grains and PAHs 
(c.f. Maloney et al. 1 19961 Bakes & Tielens ll994 . 
ii. Emission from meta-stable lines of [CI] 9823, 9850 A and 
[OI] 6300 A; fine-structure line cooling of [CII] 158 and 
[OI] 63 and 146 (im as well as Lyman a emiss ion (c.f. 
Mal oney e t al. ll996t Tielens & Hollenbach fl985l Spaans 
et al. ll994l) . 

Hi. Molecular line cooling by warm molecular gas containing 
CO, H2, H2O and OH as well as gas-grain cooling where 
warm gas is cooled at the surfaces of lower temperature 
dust grains (c.f. Neufeld et al. ll995t Spaans & Silk 2000). 

iv. Ion-molecule reactions driven by the ionisation degree of 
~ 10~ 4 maintained by the ionisation of carbon in t he FUV 
(c.f. Black & Dalgarno ll977l van Dishoeck & Black[l986). 

v. The ionisation balance of atomic gas under the influence 
of photo-ionisation reactions driven by X-ray photons and 
charge transfer. Recombination of ions on grain surfaces 
is a major ionic loss route at electron fractions les s than 
10~ 3 (c.f. Lep p & D algarno[l988; Bakes & Tielens ll994l 
Maloney et al. ll996l) . 

The global structure of any XDR is defined by a few key param- 
eters, the density nn and the energy deposition rate Hx (see 
Appendix|E) per hydrogen atom. Because the heating in XDRs 
is driven by photo-ionisation, the heating efficiency is close to 
unity as opposed to that in PDRs where the phot o-electric heat- 
ing e f ficiency is of the order of 0.3 — 1.0% ( Mal onev et alJ 
1996] lBakes &Tieleiislll994l) . Unlike PDRs, XDRs are ex- 
posed to X-rays as well as FUV photons. 
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As one moves into the XDR, X-ray photons are attenuated 
due to atomic electronic absorptions. The lowest energy pho- 
tons are attenuated strongest, which leads to a dependence of 
the X-ray heating and ionisation rates at a given point on the 
slope of the X-ray spectrum. We assume, for energies between 
0.1 and 10 keV, that the primary ionisation rate of hydrogen is 
negligible compared to the secondary ionisation rate and that 
Auger electrons contribute an energy that is equal to the photo- 
ionisation threshold energy dVoill99lh . 

The treatment is described in the appendics and follows , 
in part, the unpublished and li ttle known work bvlYanl 
Also, we extend the work of Malon ev et alJ i 19961) in terms 
of depth dependence, H2 excitation and extent of the chemical 
network. 



5. XDR test models 

In this section, we consider four models with the same energy 
inputs and densities as the PDRs in Tabled The spectral en- 
ergy distribution is of the form exp(— E/l keV). The energy is 
emitted between 1 and 10 keV and Fpuv should be replaced 
by Fx in Table This spectral shape and spectral range can 
be changed depending on the application. We take the param- 
eters for the 1 keV electron to determine the electron energy 
deposition, since these parameters do not change for higher 
energies. When the spectral energy distribution is shifted to- 
wards higher energies, the X-rays will dominate a larger vol- 
ume, since the absorption cross sections are smaller for higher 
energies. Hx l n is the most important parameter for the chem- 
ical and thermal balance, where Hx is the energy deposition 
rate per hydrogen nucleus. The abundances used are given in 
Table|3 The elements H, He, C, N, O, Si, S, CI and Fe are in- 
cluded in the chemical network. The other elements listed are 
only used to calculate the photoelectric absorption cross sec- 
tion, cr pa . 

5. 1. Heating and cooling 

In Fig.^ the different heating sources are shown as a function 
of the total hydrogen column density, Ah- All heating is done 
by X-rays, but the way it is transfered to the gas depends on the 
ionisation fraction. When the gas is highly ionised, x e ~ 0.1, 
most (~ 70%) of the kinetic energy of the non-thermal elec- 
trons goes into Coulomb heating, which is the case in Models 
1, 2 and 4 where Hx/n is high to Ah > 10 23 cm~ 2 . For 
smaller ionisation fractions, x c ~ 10~ 4 , ionisation heating as 
discussed in Sect. IB.2l is important or even dominant. In Model 
3, ionisation heating and Coulomb heating are equally impor- 
tant at Ah < 10 218 cm~ 2 . In all models ionisation heating 
dominates especially at high column densities. When the ex- 
citation of H2 is dominated by non-thermal processes, colli- 
sional quenching of H2 can heat the gas. Naively, one would 
expect this dominance to occur where most of the X-rays are 
absorbed, but for high energy deposition rates Hx /n, the tem- 
perature is high and thermal collisions dominate the population 
of the vibrational levels. Non-thermal excitation is dominant at 
low temperature, i.e., low Hx/n. 



In Fig. |2j the important cooling processes are shown as a 
function of total hydrogen column density, Ah- At high tem- 
peratures (see Fig. OJi, cooling by [CI] 9823, 9850 A and [OI] 
6300 A metastable lines dominates, as is the case in the models 
with high radiation fields, Models 2 and 4. At lower tempera- 
tures, most of the cooling is provided by the fine-structure line 
[OI] 63/im (90%), e.g., at the edge in the low-radiation field 
Models 1 and 3. In each model, gas-grain cooling dominates 
for low Hx /n. In addition, specific cooling processes can be 
important in special cases. H2 vibrational cooling dominates at 
large depths in Model 2, but in Models 1, 3 and 4 it contributes 
no more than 10%. H2 vibrational cooling is split into a radia- 
tive and a collisional part. When the excitation of H2 is domi- 
nated by non-thermal electrons, the gas is heated by collisional 
de-excitation of H2. 

5.2. Thermal and chemical structure 

In Fig.|3| we show the temperature as a function of total hydro- 
gen column density, Ah- Variations in radiation field strength 
most strongly affect the high-density models. The temperature 
at the edge differs a factor of 30 in the high-density case. Since 
X-rays penetrate much deeper into a cloud than FUV photons, 
high temperatures are maintained to much greater depths into 
the clouds. Hx /n is very important in determining the thermal 
balance. When Hx /n is larger, this results in a higher tempera- 
ture. Therefore, Model 2 has the largest temperature throughout 
the cloud. Density turns out to be important as well. Note that 
models 1 and 4 have similar incident Hx / n and therefore have 
about the same temperature throughout the cloud. 

In Figs.[5]and|4l we show the fractional abundances of se- 
lected species. Hx/n is not only important in the thermal bal- 
ance, but also in the chemistry. Therefore Models 1 and 4 with 
about the same incident Hx/n, show similar abundances. The 
most striking difference with the PDR models is that there is no 
longer a well-defined transition layer C + — > C — > CO present. 
On the contrary, both C and C + are present throughout most of 
the cloud having fractional abundances of ~ 10~ 5 — 10~ 4 . Only 
at very low Hx / n, which results in a low temperature, there is 
a partial transition to CO. The transition from atomic to molec- 
ular hydrogen is much more gradual than in the PDR models. 
A considerable amount of OH is present in all models at all 
column densities. The temperature determined by Hx /n is im- 
portant. In Model 3, OH has the largest abundance (> 10~ 6 ) 
at all column densities. In other models such large fractions are 
seen only at very high depths into the cloud. The formation of 
CO and H2O is most efficient at high densities and low Hx /n. 
Therefore, these species have large abundances throughout the 
high-density, low -radiation field Model 3. In Model 4, where 
the radiation field is somewhat higher, CO and H2O reach 
large abundances only at high Ah. At low densities, they are 
only formed at large depths into the cloud (Model 1 and 2). 
Secondary ionisations are most important for the production of 
H + . Recombination is slower at lower densities. Therefore, the 
H + fractional abundance is highest in Model 4. HCN, HCO+, 
HNC, C2H, CS and SiO have much larger abundances at high 
temperatures than in the PDR models. 
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log(N H ) [cm- 2 ] log(N H ) [cm- 2 ] 

Fig. 1. Important heating processes for Model 1 (top left), 2 (top right), 3 (bottom left) and 4 (bottom right). 



o 
o 



-16 


I Model 1 










-18 


- F x = 1 .6 

7n H =10 3 
























Gas-groin 




^Jota 


cooling 




-22 






Ol 63/irn" "-. 








_ CI 9823, 9850 






x ^\ 
















-24 






\ 

yib. col 


Cll 158fjm\\ 




-26 





6300^ s 


i 


1 




21 


22 

log(N H ) 


23 
[cm- 2 ] 


24 




-15 












-16 












-17 




Total cooling 








-18 
























-19 


= Cll 158(im 


\ 

\ 


H 2 Vib. rod'" 






-20 






\ 






-21 

-22 


; Model 3 
rF x =1 .6 
; n H =10 5 ' 5 




^01 63um 

, \ 

\ \ 
, , \ v , 






21 


22 

log(N H ) 


23 
[cm" 2 ] 


24 





-16 



o 
o 
o 



-24 

-26 



-_ _CI^9823, 985( 
Ol 6300 


- ,.7 7~r^^--^£td coo 


ing 




01 63/^m 








Gas-grain 








H 2 Vib. rod 


Cll 158/inT' 


\ , ""^^ 

\ — 


\V\ _, 


LModel 2 






\t 
'i 


- F x =1 60 






I n H =10 3 






\ 


21 


22 23 
log(N H ) [cm- 2 ] 


24 





6 



Total cooling 



E H 2 Vib, rad 



E_CII 158/im 



-20 r 

E Model 4 
-21 rF x =160 

-22 P. ", ' . i 
21 




22 23 
bqfNH) [cm 



Fig. 2. Important cooling processes for Model 1 (top left), 2 (top right), 3 (bottom left) and 4 (bottom right). 
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Fig. 3. Fractional abundances and temperature for Model 1 (top left), 2 (top right), 3 (bottom left) and 4 (bottom right). 
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Fig. 4. Fractional abundances for model 1 (top left), 2 (top right), 3 (bottom left) and 4 (bottom right). 
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6. Conclusion 

We conclude this paper by a direct comparison between the 
PDR and XDR models. To emphasise that XDRs penetrate 
much deeper into cloud volumes than PDRs, we use the same 
scale for all models. Then, it is also possible to distinguish 
between gradients in abundance, cumulative intensity, column 
density and column density ratios. XDR Model 3 is only plot- 
ted to TYh ~ 10 23 5 cm" 2 , since Hx / n becomes too small and 
no reliable results are obtained at higher column densities. 

In Fig. |5] we show for Model 4 the abundances of selected 
species. At the edge, both neutral and ionised species are more 
abundant in the XDR models, and the relative abundances also 
differ with respect to one another. In the XDR for example, the 
neutral species CH and CH2 are more abundant than CH + and 
CH^~, respectively. In the PDR, this is the other way around. 
CN and CN + are almost equally abundant at the edge in the 
PDR, while CN exceeds CN + by three orders of magnitude in 
the XDR. Although the amounts of CS+ and HCS+ are larger 
than those for CS and HCS, respectively, at the edge of the 
cloud in the XDR, the abundance difference is less than in the 
PDR. The abundance of He + is five orders of magnitude larger 
in the XDR, due to secondary ionisations. H~ is enhanced by 
three orders of magnitude, due to the higher ionisation degree. 
It is also easily seen that in PDRs the fractional abundances 
vary over many orders of magnitude, while the abundances 
in XDR Model 4 stay almost constant to a column density of 
Ah w 10 22 cm -2 , where the transition from H to H2 starts. 

In Figs. [6] and Q we show cumulative line intensities for 
fine-structure lines at every column density, i.e., the emergent 
intensity arising from the edge of the cloud to column density 
A H = ^z: 

i{z)=±- r k{ Z ')dz'. CD 

27T Jo 

Although the total [CII] 158 fim line intensity is higher in the 
XDR, the flux originating from the edge to Ah ~ 10 22 cm" 2 
is higher in the PDR except when the XDR is characterised 
by very high Hx/n values which is the case in Model 2. In 
all PDR models, all carbon is in C + at the edge, while a large 
part of the carbon is neutral in XDR Models 1,3 and 4. In all 
models, oxygen is mostly in atomic form. The [OI] 63 /im line 
intensity to Ah ~ 10 22 cm" 2 is larger in the low-density XDR 
models, which is possible due to higher electron abundances. 
The intensity is lower in the low radiation, high density XDR 
Model 3, since the temperature is higher in the PDR. For Model 
4 they are about the same, since the density where the line gets 
thermalised is almost reached. In the XDR, all line intensities 
increase more or less steadily with increasing column density. 
PDRs, however, primarily affect cloud surfaces causing more 
sudden changes. The line intensities of [CI] 609 /im and 369 
/im arise from a more or less well defined part part of the cloud 
and start to increase at column densities Ah > 10 215 cm" 2 . 
The line intensities of [CII] 158 /im are larger than those of 
[Sill] 35 /im in the PDRs except in Model 4. This is in contrast 
to the XDR models, where the [Sill] 35 /im line intensity is 
always stronger. The fact that [Sill] 35 /im lines are quite strong 



in XDRs was already noted by Malon ev et alJ ( fl993) . The line 
intensities for [Fell] 26 /im and 35 /im are larger for the XDR 
models except again for Model 3. 

In Figs. [8] and |9] we show cumulative column densities 
for selected species. They illustrate again that XDRs affect 
whole cloud volumes and PDRs create layered structures. In 
PDRs, the increase in column densities are very sudden for all 
species. For example, C and CO show this due to the very dis- 
tinct C + /C/CO transition. In the XDRs, however, the increases 
in column density are much more gradual. The only sudden 
change in XDRs is where the H/H2 transition occurs. 

In Fig.[H)] the cumulative column density ratios for CO/H2, 
CO/C, HNC/HCN, and HCO+/HCN are shown as a function 
of total hydrogen column density. The ratios for the XDRs 
are almost constant upto Ah ~ 10 22 cm" 2 , unlike those 
in PDR models. In PDRs, CO/C ratios increase by approxi- 
mately four orders of magnitude from the edge (< 10" 4 ) to 
Ah = 10 22 - 3 cm" 2 (> 1). In XDRs, this ratio is constant to 
Ah ~ 10 22 cm" 2 and then increases slowly. For each cloud 
size, while keeping the energy input the same, CO/C ratios in- 
crease at higher densities. The ratios go down for higher ra- 
diation fields. For the same density and energy input, CO/C 
is lower when the cloud is irradiated by X-ray photons, with 
the exception for Model 3 where this is only valid at Ah > 
10 217 cm" 2 . CO/H2 is somewhat more complex. When only 
the energy input is increased in PDRs, this ratio is higher when 
A H < 10 21 cm" 2 . For A H = 10 22 3 cm" 2 , the ratios are about 
the same. There is also a minimum where the H/H2 transition 
occurs. This minimum is more prominent for higher radiation 
fields. In XDRs, the CO/H2 ratio is lower when the radiation 
field is higher. In PDRs and XDRs, the CO/H2 ratios are higher 
when the density is increased. When the cloud is irradiated by 
X-ray photons, CO/H2 ratios are lower, with the exception for 
Model 3 again at A H < 10 21 5 cm" 2 . In PDR Models 1, 2, 
and 3, significant column densities for HCN, HNC and HCO + 
are reached between Ah = 10 215 and 10 22 cm" 2 . Therefore, 
the HNC/HCN and HCO+/HCN ratios discussed are for col- 
umn densities A H > 10 22 cm" 2 . In PDRs, HNC/HCN is lower 
when the density is higher. No significant changes are seen for 
different radiation fields at these columns. HNC/HCN is gener- 
ally lower for high Hx /n in XDRs. At high column densities, 
where Hx /n is low, HCN/HNC ratios are equal or somewhat 
higher than those for the PDR. HCO+/HCN and HNC/HCN 
are of the same order in PDRs, but in XDRs HCO+/HCN is 
higher in most cases. 
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Appendix A: Heating processes in PDRs 

A.1. Photo-electric emission 

In PDRs the photo-electric emission from (small) dust grains 
and PAHs is the dominant heating sou rce. W e use the analytical 
expression given by Ba kes & Tielens! Jl994l) which is given by 



10 e Gg.dust n R erg cm 3 s" 



(A.l) 



wher e Gn j„„, is the radiation field attenuated by dust absorp- 
tion jBlack & Dalgarnoll 19771) and the heating efficiency e is 
given by 

4.87- 10~ 2 



[l + 4.10-3(G T k 1/2 /nc)°- 73 ] 

3.65 • 10- 2 (T k /10 4 ) ' 7 
[1 + 2 ■ 10- 4 (G T k 1/2 /n e )] ' 



(A.2) 
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1 /2 

Note that the efficiency depends on the ratio Go T k /n e , 
which is the ratio of the ionisation and recombination rates. Ay 
is the visual extinction at optica l wavelengths caused by inter- 
stellar dust. iBohlin et alJ lll978l) relate the total column den- 

N(R) 



sity of hydrogen, A/jj 
E(B - V): 



2iV(H2) to colour excess, 



5.8 x 10 21 cm 2 mag" 



(A.3) 



E{B - V") 

The visual extinction then follows consequently: Ay = 
Z. IE(B-V) and Ay = 5.34 x l(T 22 iVH. Note that the results 
of Savag e et alJ i 19771) are often used, but in this paper only H2 
(and not H I ) is taken into account. 

A.2. Carbon ionisation heating 

At the edge of the cloud, most of the carbon is singly ionised. 
The photo-electron energy released in an ionisation is AEc = 
1.06 eV. The ionisation rate, at a certain point in the cloud, is 



given by K ion = 1.76G£ 



0, carbon ' 



. The heating rate due to the 



ionisation of carbon is then given by 



n(C) AE C 



(A.4) 



After substitution of numerical values we get the following 
heating rate for the local radiation field G' carbon : 



r c = 2.79 x 1CT 22 n(C) G 0> 



carbon 



erg cm 3 s 1 , 



(A.5) 



is the radiation field attenuated 



where this time G' n „„ 

by dust absorption (Black & Dalgarno Il977l). carbon self- 
absorption JWernerll 970) and H? Jde Jong et alfl 980). 

A.3. H 2 photo-dissociation heating 

Absorption of Lyman- Werner band photons leads to the exci- 
tation of H2. About 10% of the excitations l eads to decay into 
the continuum of the grou nd electronic state jField et all 19661 
IStecher & Williamsll 19671) . The heating related to this dissoci- 
ation is given by 



T = 0.1 Koxc. < E diss >, 



(A.6) 



where < i? d i SS > i s the mean ki netic energy of the H atoms 
and is set to 0.4 eV JSpaansll9 96). The excitation rate of H9 is 



given by k cxc . = 3.4 x 10" 
radiation field given by 



-10 r" 

U 0,H 2 



Go,h 2 =P(r)G ex P (-2.5 A 



where G H is the local 



(A.7) 



Self-shielding is explicitly taken into account for the excita- 
tion of H2, by the introduction of the shielding factor /3(t) (see 
ID.2.3> . After substitution of numerical values we get a heating 
rate of 



T H2 - 2.2 x 10~ 23 P(t) Go 

exp(— 2.5 Ay) erg cm -3 s . 



(A.8) 



A.4. H 2 collisional de-excitation heating 

FUV excitation is followed by decay to ro-vibrational levels in 
the ground state. Collisional de-excitation leads to gas heating. 
This cascade process is very complicated, but we simplify this 
process by using a two-level approximation (see ID.2.21 . The 
resulting heating rate is given by 



r H2 - KH) 7l H +n(H2) 7l H o 2 ]n(H 2 V)£;, erg cm - 3 s~Ha.9) 



where the coefficients are given by iHollenbach & McKed 
ll 19791) 



7l H = HT 12 T k a5 exp(-1000/T k ) (A.10) 
7l H 2 = 1.4 x 10" 12 T k °- 5 exp(-18100/(T k + 1200)). (A.ll) 

Both of the above expression are in units of cm 3 s _1 . 

A.5. Gas-grain collisional heating 

When gas and grains differ in temperature they can transfer 
heat through collisions. The heating rate of the gas is given by 
llHollenbach & McKeell979l fl989) 



1.2 x 10" 31 n 2 



( T k 



1/2 



100 A 



V 1000 J \ a r _ 
x[l-0.8exp(-75/T k )](T d -T k ). 



1/2 



(A. 12) 



The minimum grain size is set at a m ; n = 10 A and the dust 
temperature Td is given by 

T d = (8.9 x 10" 1 V o G o exp(-1.8A y ) + (A.13) 
2.7 5 + 3.4 x 10" 2 [0.42 - ln(3.5 x 10" 2 ri 00 T )] 
xtiooTq) 02 , 

based on the results of IHollen bach et aLN 19911) 

A.6. Gas-grain viscous heating 

Radiation pressure accelerates grains relative to the gas and the 
resulting drag contributes viscous heating to the gas. Grain ac- 
celeration time scales are short compared to other time scales, 
and therefore the grains may be considered moving at their 
local drift velocity, v^. All the momentum is transferred to 
the gas, predominantly by Coulo mb forces. For drift velocities 
Vd < 10 3 cm s" 1 (Spit zerll978h . no significant gas-grain sep- 
aration takes place. In the following we take Vd = 10 2 cm s . 
The heating rate is given by 



87re 4 n d ^d( fc7 k) _1 (lnA)w d 
[n(C+)G(y c+ )+n e G(y e )}, 



(A. 14) 



where n d is the grain volume density, Z d is the grain charge, 
n(C + ) and n c are the respective C + and electron volume den- 
sities and the functions A and G(y) are given by 
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A 



1.5 Z A 



\kT k )^(7rn e ) 



-0.5 



G(y) = ^{erf(y) 



7vjy e y }» 



(A.15) 
(A. 16) 



where y — Vd/vth and u tn the thermal velocity of C + ions and 
electrons. The error function erf (y) is given by 



crf (y) = / ' 
lo 



dt 



(A. 17) 



A. 7. Cosmic-ray heating 

At large column densiti es, co s mic ra y he ating can become im- 
portan t. IGlassgold & Lang erl dl973h and ICravens & Dalgarriol 
(1978) calculated that the amount of heat deposited in a 
molecular gas is about 8 ke V per primary ionisation. Then, 
iTielens &*H ollenbach ill 9851) find for the total heating rate, in- 
cluding helium ionisation 



r CR = 1.5 x 1CT 11 C n(H 2 ) erg cm" 3 s _1 , 



(A.18) 



where £ is the cosmic ray ionisation rate per H2 molecule. 

Appendix B: Heating processes in XDRs 

B. 1. Heating due to Coulomb interactions 

When X-rays are absorbed, fast electrons are produced. These 
fast electrons lose part of their energy through Coulomb inter- 
actions with thermal electrons, so the X-ray heating is given 
by 



B.2. Heating due to H 2 ionisation 

H? i onisation can lead to gas heating ( IGlassgold & Langerl 
1973). When H2 is ionised by a fast electron and subsequently 
recombines dissociatively, about 10.9 eV (1.75 x 10~ n erg) of 
the ionisation energy can go into kinetic energy. H^ can also 
charge transfer with H. This is an exothermic reaction, with 
an energy yield of 1.88 eV, of which we assume half, 0.94 eV 
(1.51 x 10~ 12 erg), to go into heating. H^ can also react to H^", 
and subseque ntly recombine dissoc iatively or react with other 
species. Glassgold & Langerl Jl973l) argued that for every Hj 
ion formed 8.6 eV (1.37 x 10 -11 erg) goes into gas heating. 
The H2 ionisation rate cooling is then given by 



17.5fc r 



H2 ion 



1.51fcH£H + 13.7fcn 2 a;H2 



10 -12 £h 2 x h 2 n er S cm _3 s _1 , 



(B.5) 



where k c , fcn and fcjj 2 are the rates of dissociative recombi- 
nation, charge transfer with hydrogen and the reaction to H^", 
respectively. 

B.3. Gas-grain collisional heating 



We use the results of Sect. IA. 51 The dust temperature was found 
bv lYar1lll997l) : 



T d = 1.5 x 10\H x /x d ) - 2 K, 



(B.6) 



where x d = 1-6 x 10~ 
erg s -1 . 



is the grain abundance and Hx in 



T x = V n H x , 



(B.i) B.4. H 2 vibrational heating/cooling 



where 77 is the heating efficiency, depending on the H2/H 
ratio and the electro n abundance x. We use the results of 
Dalga rnoeta Tl lll999h . Their calculated heating efficiency 77 in 
an ionised gas mixture is given by 



V 



10r77H 2 Hc + ??HcH 

10r + 1 



(B.2) 



where r = n(H2)/n(H). 77h 2 h and t^hcH are the heating effi- 
ciencies for the ionised pure He and H2 mixture and the He and 
H mixture, respectively. Both are parametrised through 



When the vibrational levels of H2 are populated by non-thermal 
processes, thermal collisional quenching and excitation can re- 
sult in a net heating despite downward radiations. When non- 
thermal reactions are not important, H2 can be an important 
coolant. The resulting collisional vibrational heating or cool- 
ing is given by 



TH 2 vib,col — ^~>vjTlvj X 



vyC{vj 



(B.7) 



S 



vf) x 



E v iji) erg cm 3 s 



r/ = l + (77o-l)/(l + cO. 



(B.3) 



The va lues of rj , c and a are given in Table 7 of Dal garno et alJ 
Jl999l) . and x is the electron fractional abundance. It has to be 
modified when the H2-He mixture is considered: 



Where C(vj — > v'j') is the total collision rate from level vj to 
v'f in units of s . Radiative cooling due to downward decay 
of the vibrational levels is given by 



-^H 2 vib,rad 



T, V jA(vj — > v'j')n V j erg cm 3 s 1 



(B.8) 



1.83a; 
1 + 0.83x' 



(B.4) 



The population of the vibrational levels is discussed in Sect. 

DH3 
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Appendix C: Cooling processes 

C. 1. Fine-structure line cooling 

Since most of the gas is atomic in the radical region, the 
dominant coolants are the atomic fine-structure lines. The most 
prominent cooling lines are the [CII] 158 /im and [OI] 63 fim 
and 146 /im lines. For the calculation of the thermal balance 
we also take into account Si + , C, Si, S, Fe and Fe + . We use a 
compi lati on for the collisional dat a from Sternberg & Dal garnol 
|1995|) . iHollenbach & McKed dl989h. ISampson et all 
i 1994. iDufton & Kingston! i 1994. Ijohnson et alJ ll 19871) 
iRoueff & Le Bourlod Jl99oT. ISchroder et alJ dl99lh iMendozal 
Jl983h . IChambaud et all Jl980l) and Ijaauet et alJ Jl992h We 
take into account collisions with electrons, H + , H and H 2 
(ortho and para) for the excitation of the species to different 
levels. In the PDRs, collisions with H + are not the dominant 
excitation source but in XDRs the ionised fraction of hydrogen 
can be as large as ten percent and become important for the 
excitation of some levels. 



C.2. Metastable-line cooling 



We included the 
Si, Si+, O, Q+, 
taken 



metastable 
S, S+ Fe 



from 



cooling lines of C, C , 
and Fe + . All the data is 
except for Si + 



'1989) 



Hollenbach & McKee ( 
(Dufton & Kingston 19941 C+ JSampson et all 1994 and 0+ 



McLaughlin & Bell 199 



1' 



C.3. Recombination cooling 

At temperatures higher than ~ 5000 K, cooling due to recombi- 
nation of electrons with grains (PAHs) is important. The cool- 
ing depends on the recombination rate which is proportional to 
the product n c rtn- The cooling rate increases when GqT® 5 /n e 
goes up, due to an increase in cha rge and hence Coulomb in- 
teraction. |Bake^^^elens| |l994) calculated numerically the 
recombination cooling for a variety of physical conditions. An 
analytical fit to the data is given by 



A = 3.49 ■ 10" 30 T k Q (G T k 1/2 /no) /3 (CI) 
n e nn erg s _1 cm -3 

where a = 0.944 and = 0.735/T k 068 . 

C.4. Molecular cooling by H 2 , CO and H 2 O 

For the rotational and vibrational cool ing of H 2 , CO and 
H 2 0, we use the fitted rate coefficients o flNeufeld & Kaufmanl 
( 1993b and lNeufeldetail l l 19951) . They present a cooling rate 
for species i through: 



-3 ,,-1 



A = L n(xi) n(H 2 ) erg cm s 



(C.2) 



where n(H 2 ) and n{xi) are the densities of H2 and species Xi, 
respectively. L is given by 



1 

I 



1 



n(H 2 ) | 1 
Llte Lo 



n(H 2 ) 



ni/2 



1 



H1/2 



'LTE 



(C3) 



We in terpo late in the tables g iven by iNeufeld & Kaufmanl 
(1993) and Neufel d et alJ |l995), to find the values L , n 1/2 
and Llte and a. Lo is the cooling rate coefficient in the low 
density limit and n 1 / 2 is the H2 density where L has fallen by 
a factor of two below Lq- a is chosen to minimize the maximal 
fractional error in the fit at other densities. Lq is a function of 
temperature, and Llte, ^i/2> an d a are functions of temper- 
ature and the optical depth parameter N(xi), which is given 
by the gradient N(xi)/5vd- N(xi) is the column density of the 
species 2^. To take into account collisional excitation by elec- 
trons and atomic hydrogen, we follow Yan ( 1997) and replace 
n(H 2 ) by n Iot and n V it>- F° r H 2 rotational and vibrational cool- 
ing, n rot and n v ;b are given by 



ra ro t(H 2 ) = n vib (H 2 ) = n(H 2 ) + 7n(H) + 16n(e). (C.4) 
For rotational cooling by CO, n lot is given by 



n rot (CO) - n(H 2 ) + 1.414n(H)a H /(TH 2 
+1.3 x 10~Vc)/cr H2 -u, 



(C5) 



where cr H = 2.3 x 10~ 15 crrT 2 , cr H2 = 3.3 x 
lO-^Tk/lO 3 ))- 1 / 4 cm" 2 and v = 1.03 x lO 4 !^ 5 cms" 1 . 
For H 2 rotational cooling, n rot is given by 



n rot (H 2 0) = n(H 2 ) + 10n(H) + 

n(e)fc e (l,20,1.9,T k )/fc H2 , 



(C6) 



where k H . 2 = 7.4 x lO^T^cm^- 1 and k c {i,b,d,T k ) are 
the H 2 and electron impact excitation rate coefficients, respec- 
tively. k e (i, b, d, T k ) for the excitation from level i — > i + 1 in 
units of cm 3 s _1 is given by 



k c (i,b, d,T k ) 



3.56 x 10~ 6 d 2 
T k °- 5 [2-l/( l + l)] 



In 



CAE 



C 



exp(f3AE) x (C.7) 



-0.577 \ 



6XP [l + 2(3AEJ 



where b is the rotational constant in cm 1 , d the dipole moment 
in Debye, (3 = 11600/T k , AE = 2.48 x 10" 4 fe(i + 1) and C 
is given by 



C 



c = 



9.08 x 10 3 

b(i + l) 
1.93 x 10 4 



d < 1.53 



(C8) 



— exp(-1.18/d 3 ) d>1.53 



db(i + 1) 

For CO vibrational cooling, n v ;b is given by 
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n vib (CO) = n(H 2 ) + 50n(R) + 7i(e)L co ,e/ X co ,o (C.9) 
where 

L COe = 1.03 x 10" 10 (T k /300) a938 exp(-3080/T k )(C.10) 
L C o,o = 1.14 x 10- 14 cxp(-68.0/T k 1/3 )exp(-3080/T k ) 

For H 2 vibrational cooling, n v ib is given by 

n vi b(H 2 0) = n(H 2 ) + 10n(H) + n(e)L H2 o,e/L H2 o,o (C.ll) 
where 

Lh 2 o,c = 2.6 x 10" 6 T 1 7°' 5 exp(-2325/T k ) (C.12) 
Lh 2 om = 0.64 x 10" 14 cxp(-47.5/T k 1/3 )exp(-2325/T k ) 

In the XDR models, H 2 vibrational cooling is treated dif- 
ferently, since non-thermal processes play an important role, 
which is discussed in Sect. IB. 4l 

C.5. Cooling by electron impact with H 

The cooling due to the excitation of hydrogen is imp ortant at 
tempe ratures T > 5000 K. The cooling rate is given bv lSpitzeil 
lll978l) : 



i?H 2 = 0.5 n H «h n d er d e R2 S R (T k ) (D.l) 
« 6 x 10~ 17 (T k /300) a5 n H n en 2 S(T k ) cm" 3 s" 1 , 

where and Od are the volume density and cross section 
of dust grains and njj, i>h and S(T k ) are the volume density, 
thermal velocity and thermally averaged sticking coefficient 
of hydrogen atoms. We use the sticking coefficient given by 
iHollenbach & McKed J 19791) 

S(T k ) = [1 + 0.4(T k + T d )°- 5 + 2 x 10~ 3 T k + (D.2) 

8 x 10~ 6 Tg]- 1 , 

where Th is the dust temperature. E q. JD.2t is the same as eq. 
(4) in Tielens & Hollenbach ( 1985), except for the term ejj 2 , 
the recombination efficiency, which is given by 

V 2/fe 2 otpc J 

where fi is the H 2 fraction that stays on the surface after for- 
mation, /?h 2 an d Pa are the desorption rates of molecular hy- 
drogen and physisorbed hydrogen atoms, respectively, F is the 
flux of hydrogen atoms and a pc is the evaporation rate from 
physisorbed to chemisorbed sites. The se three terms do minate 
in dif ferent temperature regimes. See lCazaux & Tielensl J20021 
2004) for a more detailed discussion. 



A 6 _h = 7.3 x 10~ 19 n c n(H) (C.13) 
x exp(-118400/T k ) erg cm _3 s -1 . 

Appendix D: Chemistry 

For most of the chemical reaction rate s, we make use of th e 
UMIST database for astrochemistry bv lLe Teuff et all ( 2000). 
In de PDR model we use a network containing all the species 
with a size up to 6 atoms. For the XDR model we use all species 
with sizes up t o 3 atoms an d some of 4 atoms. These species 
are taken from lYanUl997l) . Below we discuss the additional 
reactions. 

D.1. H 2 formation on dust grains 

The formation of H 2 is very efficient o ver a wide r ange 
of tem peratures. It was already shown by iGould & Saroeter 
(1963) that H 2 is not formed efficiently in the gas phase. 
Most of the formation, which is still n ot very well understood , 
takes pla ce on grain surfaces JHol lenbach & Salpete3ll97ll) . 
Recently, ICazaux & Tielensl J2002L 120041) developed a model 
for the formation of hydrogen under astrophysically relevant 
conditions. Th ey compared their resu lts with the laborator y 
experiments bv IPirronello et all |l999) and lKatz et all |l999). 
They find a recombination rate of 



D.2. Additional reactions in PDR models 
D.2.1. Recombination on PAHs 

Collisions of electrons and ions with grains can become an 
important recombination process in dense clouds of low ion- 
isation. We in clude reactions with PAHs following Sect. 5 of 
IWolfire et all J2003I) in the PDR models. The C+/C transition 
occurs at larger column densities when PAHs are included. 
Deep into the cloud, the electron abundance is reduced by sev- 
eral orders of magnitude. 

D.2. 2. Vibrational excitation of H 2 

In PDRs, molecular hydrogen can be excited by absorption of 
FUV photons in the Lyman- Werner bands. Fluorescence leads 
to di ssociation in about 10 % in of the cases (see Field et al. 
119661 Stecher & Williams 1 19671) . and in the remaining 90% 
of the cases to a vib rationall y excited state of the ground elec- 
tronic state (Black & Dalg arnoll976l) . To simplify matters, we 
treat the electronic ground state as having a vibrational ground 
state and a single excited vibrational state. iLondonl dl978l) 
found that the effective quantum number for this pseudo-level 
is v = 6, and the effective energy is E*/k = 2.6 cV/k = 
30163 K. We treat excited molecular hydrogen, H 2 V, as a sep- 
arate species in our chemistry. H 2 V can be destroyed by direct 
FUV dissociation, radiative decay or collisional de-excitation, 
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and chemical reactions with other species. Since vibrational de- 
cay is a forbidden process, a large abundance of H2V can be 
maintained. H 2 V can react with other species with no activa- 
tion barrier or a reduced one. In the UMIST database, the rates 
for a reaction between two species are parameterised as 

R = a (T k /300)' 3 cxp(-7 /T k ) cm 3 s" 1 . (D.4) 

For reactions with H2V, 7 is replaced by 7* = max(0.0, 7 - 
30163). When reactions have an activation barrier lower than 
2.6 eV, the barrier is set to zero. When t he barrier is la rger than 
2.6 eV , the barrier is reduced by 2.6 eV. lTielens & Hollenbach 
( 1985) state that for important reactions such as 

H 2 V + C+ -> CH+ + H 
and 



We assume that Si and CI get singly ionised after absorbing an 
X-ray photon. We also include C 2+ , N 2 +, S 2+ and Fe 2 +. The 
direct (or primary) ionisation rate of species i at a certain depth 
z into the cloud is given by 

Ci, P rim = / ai(E) F ^ Z ^ dE, (D.5) 

JB miB E 

where the ionisat i on cross sections cr, are taken from 
IVerner & YakovlevNl995l) . 

D.3.2. Secondary ionisations 

Part of the kinetic energy of fast photoelectrons is lost by ioni- 
sations. These secondary ionisations are far more importan t for 
H, H2 and He than direct ionisation. Dalgarno etaDCi^ cal- 
culate the number of ions A^ on produced for a given species i. 
For a given electron energy E, N- lon is given by 



H 2 V + O -> OH + H, 

this is a good approximation since the activation barrier of ~ 
0.5 eV is a lot smaller than the vibrational excitation energy of 
2.6 eV. For reactions with barriers of the same order or larger 
one can overestimate the reaction rates. 



N ion = E/W, (D.6) 

where W is the mean energy per ion pair. iDalgarno et alJ 
( 1999) calculated W for pure ionised H-He and H 2 -He mix- 
tures and parameterised W as: 



D.2.3. Shielding of H 2 and CO 



W = W (l + cx a ), 



(D.7) 



In PDRs, the photo-dissociation rate of both H 2 and CO is in- 
fluenced by line as well as continuum absorption. The disso- 
ciation rate of H 2 is decreased by self-shielding. For an H 2 
line optical d epth r < 1 0, we adopt the self-shielding fac- 
tor given by dShulllll978h . When the line absorption is dom- 
inated by the Doppler cores or the Lorentz wi ngs (i.e., r > 
10), w e use the self-shielding factor as given bv lde Jong et all 
(1980). The CO photo-dissociation rate is decreased by both 
CO self-shielding and H 2 mu tual shielding. We use Table 5 of 
Ivan Dishoeck & Blackl Jl988l) . to determine the shielding fac- 
tor as a function of column densities iV(H 2 ) and N(CO). 

D.3. Additional reactions in XDRs 
D.3.1. Primary ionisations 

In the XDRs we do not use the photo-ionisation rates from 
UMIST. X-rays are absorbed in K-shell levels releasing an 
electron. An electron from a higher level may fill the empty 
spot and with the energy surplus another so called Auger elec- 
tron is ejected. This process leads to multiply ionised species. 
Due to charge transfer with H, H 2 and He, they are quickly 
reduced to the doubly ionised state. We therefore assume that 
the ionisation by an X-ray photon leads to a doubly ionised 
species, as does absorption of an X-ray photon by a singly 
ionised species. When rates for charge transfer with H and He 
are very fast, elements are quickly reduced to singly ionised 
atoms, which is the case for 2+ , Si 2+ and Cl 2+ . Therefore, 
we add only 2+ to the chemical network to represent them. 



where Wo, c and a are given in Table 4 of their paper. The 
corrected mean energies for ionisation in the H-H 2 -He mixture 
are given by 



W(H+) = W H ,He(H + 

W(H+)=Wh 2 ,hc(HJ) 



[1 + 1.89 
1 + 0.53 



n(H 2 ) 



n(H) 
n(H) 



n(H 2 ) 



(D.8) 
(D.9) 



The ionisation rate at depth z into the cloud for species i is then 
given by 



-Emin 



a p3 ,(E)F{E, z)—dE s" 1 per H nucipMO) 



We rewrite this to a rate dependent on the fractional abundance 
of the species xf. 



E 



, a pa (E)F(E,z)——dEs 1 per specpfeK}) 

-Emin W Xi 



where Xi is the fraction of species i. Since we integrate over 
the range 1-10 keV and W goes to a limiting value, we use the 
parameters applicable to the 1 keV electron. The ionisation rate 
then simplifies to: 
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1 keV 



W(l keV)xi 

1 keV 
W(l keV)x. 



a pab (E)F(E,z)dE (D.12) 



■Hx s per species i. 



We also include secondary ionisations for C, N, O, Si, S, CI, 
Fe, C + , N + , + , S + and Fe + . We scale the ionisation rate of 
these species to that of atomic hydrogen by 



Q = Ch- 



(D.13) 



We integrate over the range 0.1-1.0 keV to get an average value 
of the electron impact i onisation cross sect ion <t c1 . Using the 
experimental data fits of lLennon et alJ dl988h . The scaling fac- 
tors (Tcm/ctclh for C, N, O, Si, S, CI, Fe, C+ N+, 0+, S+ 
andFe+ are 3.92, 3.22, 2.97, 6.67, 6.11, 6.51, 4.18, 1.06, 1.24, 
1.32, 1.97, and 2.38, respectively. 

D.3.3. FUV photons from secondary electrons 

When energetic electrons created in X-ray ionisations collide 
with atomic and molecular hydrogen, H2 Lyman- Werner and H 
Lyman a photons are produced, which can significantly affect 
the chemistry. The photoreaction rate Ri per atom or molecule 
of species i is given by 



Ri 



1 - LJ 



(D.14) 



The values of p a are taken from table 4.7 o f lYanl ill 9971) and val- 
ues of p m are the r ates for cosmic-ray induced reactions from 
Ee Teuff et al. (2000). There is an exception for CO, however, 
where we take the ra te, corrected for self-shielding, given by 
iMalonev et all J 19961) : 



yield to higher levels decreases very rapidly. First we calculate 
the mean energy for excitation, W . in the H?-He mixture . The 
parameters are listed in Table 5 of Dalgarno et al. ( 1999j). The 
function W has the same form as equation ID.7I . The mean 
energy for excitation also depends on the abundances of H and 
H2. The yield has to be corrected with a factor C(H,H2), which 
is given by 



C(H,H 2 ) = 
C(H,H 2 ) = 



2n(H 2 ) 



n(H) + 2n(H 2 ) ' 
n(H 2 )/rc(H) 
n(H 2 )/n(H)+o(x) 



x > 10" 



(D.16) 



10 



7 < x < 10" 4 



where a (x) = 0.5(x/10 4 ) 015 . The rates f or excitation by 
thermal electrons are taken from lYarld 19971) . who finds that 
the transitions rates for H 2 (w = 0) to H 2 (i> = 1, 2) are given 
by 



R(0 -> 1) = 9.7 x 10" n (T k /300) - 87 exp(-6140/Tiip.l7) 



R(0 



7.5 x 10-^(T k /300) u - yi exp(-11900/S).18) 



The excitation rate for the transition v — ► v + 1 is taken to 
be v times the — > 1 rate. Excitations with Av > 1 are not 
taken into account. The quenching rates are calculated through 
detailed balance. The quenching ra tes from v — > v' by atomic 
hydrogen are given in table 4.2 of lYanl dl997l) and are of the 
form: 



R{v 



= a(T k /300) /3 exp(- 7 /T k ) 



cm 3 s 1 . 



(D.19) 



The excitation rates are obtained by detailed balance. For the 
molecular excitation and quenching rates we use the results of 
iTine etalJdl997l) . Collisions where either before or after one 
of the H 2 molecules is in the v = state are considered. The 
rate coefficients are of the form: 



R co = 2.7x C o /2 (Tk/1000)°- 5 C H2 £H 2 s 



D.3.4. Vibrationally excited H s 



(D.15) 



Vibrationally excited H 2 can enhance reactions with an activa- 
tion barrier and also be an important heating or cooling source. 
To calculate the populations of the vibrational levels of H 2 , we 
take into account: 

- Collisions with fast electron produced by X-ray photo- 
ionisation. 

- Collisions with thermal electrons, H, H 2 and He. 

- Chemical destruction and production in chemical reactions. 

- Radiative decay. 



We use the results of Dalga rno et alJ d 19991) to calculate the X- 
ray induced excitation to the vibrational levels v = 1 and v = 
2. The ratio of the yields Y(v = 2)fY(v = 1) is about 0.070. 
Excitation to higher levels is not taken into account, since the 



log 10 R{v u v' i; v 2 ,v' 2 ) = A + B/T k + Clog 10 T k , (D.20) 



and are given in table 1 of ITine etal 1 d 19971) . who also con- 
sidered collisions with He. They give a rate coefficient for the 
v — 1 — > transition: 



log 10 i?(1^0) 



-8.8T k 
-18.9T, 



= -47.4T, 



-1/3 _ 
-1/3 
-1/3 



16.5 T k <90K(D.21) 

- 14.2 90 < T k < 230K 

- 9.4 T k > 230K. 



For the other transitions with Av = 1, the same rates are used. 
The upw ard transitions can be obtained by detailed balance. 
lYanldl997l) also calculated the dissociation and ionisation rates 
by thermal electrons and since the ionisation threshold is much 
higher than the vibrational energies one rate is used for all vi- 
brational energies: 
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R e ,di ss = 7.03 x 10~ 9 (T k /300)°- 41 exp(-118600/T k ) (D.22) 
R c ion = 8.9 x 10" 10 (T k /300)°- 57 exp(-179400/T k ) (D.23) 

The dissoci ation r ates by atomic hydrogen are given in table 
4.3 oflYar3 Jl997h . which are of the same form as equation 
( ID. 19i . F or the dis sociation rates by H2 we use the results of 
iLepp & Shulll Jl983l) . which are given by 

i?H 2 ,d iss = 6.29 x 1(T 15 x (D.24) 
Aexp(lMv - 0.037w 2 )/(T k )//(4500K), 

where A = 1.38, /(T k ) = T k 05 a exp(-a), a = [1 + (E th + 
l)/fcT k ] and E th = 4.48eV - E(v). For the dissociative at- 
tachment reaction: 



we use the results of IWadehra & Bardslevl dl978l> and the 
reaction rates have the same form as equation ( ID. 191 . 
Vibrationally excited H 2 can be destroyed in chemical reac- 
tions. Endothermic reactions with vibrationally excited H 2 can 
lower the activation barrier, by using the energy of the vibra- 
tional level. The barrier is reduced, but cannot become nega- 
tive: E' = min(0, E — E(v)). When H 2 is formed in chemical 
reactions which are exothermic, part of the formation energy 
goes into the excitation of the vibrational levels. Formation of 
H 2 on grains can play a very signific ant role. H 2 ha s a bin ding 
energy of 4.48 eV. Following Sternberg & Dalgarno ( 1989), we 
assume one third of this energy to be distributed statistically 
over all the vibrational levels: 



The interval [-Emin^max] is the spectral range where the en- 
ergy is emitted. The photoelectric absorption cross section per 
hydrogen nucleus, cr pa , is given by 

a pa (E) = J2 ACtotalViCE). (E.2) 

i 

M orrison & McCammonl dl983h state that the X-ray opac- 
ity is independent of the degree of depletion onto grains. 
Therefore, we take the total (gas and dust) elemental abun- 
dances, Ai (total), as given in Table |2] to calculate cr pa . 
The X-ray absorp t ion cr oss sections, er^, are taken from 
IVerner & Yakovlevl dl995h . The flux F(E, z) at depth z into 
the cloud is given by 

F(E, z) = F(E, z = 0) exp(- f 7 pa (S) N u ), (E.3) 

where iVn is the total column of hydrogen nuclei and F(E, z — 
0) the flux at the surface of the cloud. 



cxp(-£fo)/1.493) 
Ij v exp(— E(v)/1. 493) 



(D.25) 



where x(H 2 (v)) is the fraction of H 2 formed in vibrational state 
v. When H 2 formation reactions are endothermic, all the H 2 
is in the ground vibrational state. When they are exothermic 
part of the energy is distributed statistically following equa- 
tion ( ID.25> . The Einstein A coe fficients for radiative decay are 
taken from iTurner et alJ Jl977h . We take a weighted average 
over the rotational levels of each vibrational level, which we 
assume to be thermalised, to get an Einstein A coefficient for 
the decay from v — » v'. 

Appendix E: Energy deposition rate per 
hydrogen nucleus 

The photon energy absorbed per hydrogen nucleus, Hx, is 
given by 



p £" max 

H x = a pa (E)F(E,z)dE. (E.l) 



